Statistics of Lagrangian quantum turbulence 
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We consider the dynamics of small tracer particles in turbulent quantum fluids. The complicated 
interaction processes of vortex filaments, the quantum constraints on vorticity and the varying 
influence of both the superfluid and the normal fluid on the tracer particle effectively lead to a 
superstatistical Langevin-like model that in a certain approximation can be solved analytically. An 
analytic expression for the PDF of velocity v of the tracer particle is derived that exhibits not 
only the experimentally observed v~ 3 tails but also the correct behavior near the center of the 
distribution, in excellent agreement with experimental measurements and numerical simulations. 
Our results are universal and do not depend on details of the quantum fluid. 



Quantum turbulence is a phenomenon of utmost in- 
terest in current fluid mechanics research [l|-|lOl|. The 
turbulent behavior of a quantum liquid such as 4 He is 
very different from classical turbulence since vortices are 
quantized. This means the circulation cannot take ar- 
bitrary values as in classical turbulence, and there is 
also no viscous diffusion of vorticity as in classical tur- 
bulence. Recent measurements [l| and simulations 0, Q 
have shown that this has profound influence on various 
measurable observables, most notably the velocity distri- 
butions of small test particles embedded in the turbu- 
lent flow. Whereas for classical turbulence there is near- 
Gaussian behavior pd| - [l3j . one typically observes power 
laws for quantum turbulence [lp]!, [l4| . 

The velocity statistics has been subject of several re- 
cent papers. Paoletti et al [l| conducted a seminal ex- 
periment using solid hydrogen tracers in turbulent su- 
perfluid 4 He and found that the distribution of velocity 
components v% of the tracer particles exhibits a power 
law p(vi) oc v~ 3 distribution for large values of Vj. Simi- 
lar results were confirmed by White et al They per- 
formed numerical simulations of quantum turbulence in 
a trapped Bose-Einstein condensate by calculating solu- 
tions of the Gross-Pitaevski equation. The associated 
PDF of each velocity component Vi(i = x, y, z) was com- 
puted directly and it was confirmed that the velocity 
statistics is non-Gaussian and obeys a power-law distri- 
bution p(vi) oc v^ b with —3.6 < b < —3.3. In the fol- 
lowing, for ease of notation, we often suppress the index 
i. Adachi et al 3] numerically computed the velocity 
field of a superflow by calculating the Biot-Savart ve- 
locity induced by vortex filaments in steady countcrflow 
turbulence. They found that the resulting PDF exhibits 
a near-Gaussian distribution in the low-velocity region 
whereas a power-law p(v) oc v~ 3 is observed in the high 
velocity region. 

Apparently there is clear evidence from numerical and 
experimental approaches that power laws in the veloc- 
ity statistics are highly relevant in quantum turbulence, 
and that typically the observed power law exponent is 
close to 3. What is missing so far, however, is a theory 
of Lagrangian quantum turbulence, by which we mean a 
theory that consistently describes the dynamics of tracer 
particles of a given size embedded in the quantum tur- 



bulent flow, which would be the proper theoretical tool 
to explain the observed velocity distributions. Whereas 
Lagrangian turbulence is a well-established subject area 
for classical turbulence [l5l - [2lj |. very little is known for 
the quantum case. 

In this paper we will introduce a simple but power- 
ful dynamical model of the dynamics of a tracer parti- 
cle embedded in a quantum liquid. This model will be 
based on a superstatistical stochastic differential equa- 
tion. The superstatistics concept, introduced in [22| . has 
proved to be a very powerful method for modeling a vari- 
ety of complex systems (23l - [26j , including driven nonequi- 
librium situations and classical hydrodynamic turbulence 

Here, for the first time, we apply this concept to quan- 
tum turbulence. The result is a dynamical theory that 
quite precisely reproduces the observed velocity statistics 
in quantum turbulence and that also allows for some an- 
alytic predictions. In particular, the power law exponent 
—3 follows from our theoretical consideration in a natural 
way, and moreover a universal prediction for the entire 
shape of the velocity distribution is obtained, which is in 
excellent agreement with experimental measurements. 

Let us denote the velocity of a Lagrangian tracer par- 
ticle embedded in the quantum liquid by v(t). We start 
from a simple local dynamics which will later be extended 
to a superstatistical model. Consider a linear stochastic 
differential equation of the form 



v(t) = -Tv(t) + ZL(t) 



(1) 



Here L(t) is a rapidly fluctuating stochastic process rep- 
resenting rapid forces in the quantum liquid on a fast time 
scale, and T and £ are 3x3 matrices. The above equa- 
tion simply says that locally a tracer particle is driven 
by chaotic forces L{t) from the turbulent flow and at 
the same time there are damping processes, described by 
F. Since the chaotic forces act rapidly we approximate 
L(t) by Gaussian white noise. T and E are matrix- valued 
stochastic processes which evolve on a much larger time 
scale than L(t). The particle is driven by a mixture of 
normal and superfluid, and depending on which compo- 
nent dominates, the effective friction described by T will 
be very different. 

A characteristic property of quantum turbulence is 



2 



a spatio-temporally varying vorticity field represented 
by 1-dimensional topological vortices that reconnect and 
merge at random moments of time. A test particle may 
rotate for a short while around a local unit vector e whose 
direction will be a random variable, describing a given 
vortex filament in the quantum liquid. Hence, as a spe- 
cial case of eq. ([T]) we may consider the local dynamics 



-j(t)v +u[e(t) x v] + aL(t) 



(2) 



We assume that the damping constant 7 and the noise 
strength a are functions of t, and so is oj and the 
direction of e. The second term on the right hand side of 
eq.© represent the rotational movement of the particle 
around the vortex filament. The unit vector e and the 
noise strength a evolve stochastically on a large time 
scale T e and T a respectively. 



vortex filament is given by Q 
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(6) 



where n = ss 9.97 x 10~ 4 cm 2 /s is the quantum of 
circulation, h is Planck's constant and m is the mass of 
the fluid atom, in our case helium. 

If the tracer particle comes close to a vortex filament, 
it will typically follow a circular path around the vortex 
filament, with v = = j 1 -, where T is the period of 
one rotation. Note that the angular frequency entering 
eq. dU is thus uj = ^ = 

For an ordinary spherical Brownian particle in a 
viscous liquid one has constant damping due to Stoke's 
law: 



A special coordinate system would be e = (0,0,1), 
then ex v = (—v v , v X) 0) and the velocity components of 
the particle satisfy 



v x = ~iv x - tovy + aL x (t) 

Vy = — fUy + iOV x + <jLy(t) 

v z = -~fv z + <jL z (t) 



(3) 



If we introduce a complex variable z by defining z = 
v x + iv y , then the (x, i/)-dynamics can be written as 

z = v x + iv y = (—7 + iui)z + a(L x + iL y ) (4) 

Forming the average (• • • ) over all realizations of the noise 
L(t) one obtains on a time scale where 7 and uj are suf- 
ficiently constant 

(*(*)) = z(0)e" 7 *(cos(o;i) + i sin(wt)) (5) 

which is just damped spiraling motion around a local unit 
vector with frequency uj. We remind the reader that the 
basic idea of the superstatistics approach is to regard the 
parameters of a local stochastic differential equation as 
random variables as well [27| . This means both 7 and uj 
can take on very different values during time evolution, 
and so can the direction of e. A very small 7 corresponds 
to nearly undamped motion for a limited amount of time. 
A very small uj corresponds to almost no rotation, i.e. 
straight movement for a limited amount of time. All 
these cases are included as possible local dynamics and 
averaged over in the superstatistical approach. 

In a quantum turbulent flow, the superfluid com- 
ponent flows without dissipation while being subject 
to certain quantum mechanical constraints. These 
quantum restrictions imply that the typical form of 
rotational motion allowed in the superfluid component 
is in the form of a thin vortex line, whose circulation 
around its core is quantized rather than arbitrary as 
in classical fluids. The magnitude of the velocity field 
of the fluid particle at distance r from the core of the 
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Here v is the kinematic viscosity of the liquid, p is the 
fluid density, M is the mass and a the radius of the 
tracer particle. 

For quantum turbulence, the effective dissipation act- 
ing on the tracer particle is influenced by many com- 
peting effects, and it fluctuates strongly depending on 
whether the particle is close to a vortex filament or not. 
Far away from a vortex filament, the movement will be 
dominated by Brownian motion similar as in a normal 
liquid, whereas close to a vortex filament the movement 
will be very rapid and almost friction free, dominated by 
the superfluid. 

To take into account the fact that the effective friction 
in eq. ([2]) is fluctuating, we may write quite generally 
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(8) 



where L is a characteristic length scale and the X, are 
dimensionless random variables that evolve in time and 
space. We have squared the random variables because for 
physical reasons 7 must always be positive, though values 
close to are possible, n denotes the number of degree of 
freedoms that influence the fluctuating effective friction. 
Of course, the simplest model is to assume that the Xi 
are a rescaled sum of many microscopic random variables 
that act almost independently. Thus the Central Limit 
Theorem suggests to assume that the X, are Gaussian 
random variables. 

The quantum mechanical constraint given in Eq. ([5]) 
tells us that the average rotational velocity of the tracer 
particle is very high near the vortex core (for small dis- 
tance r). Therefore, the effective viscosity 7 acting on the 
tracer particle in eq. (J5|) is small if the particle is very near 
to the vortex core. This means that ^™ Xf is small. On 
the other hand, if the test particle is very far from the 
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vortex filament, then 7 is large and the friction effects 
are strong, mainly due to the normal fluid component. 
This suggests the physical interpretation that the Xi's 
may just be identified with the perpendicular distances 
of the test particle from the nearest vortex filament. The 
vortex filaments themselves of course evolve in a highly 
complicated stochastic way. Since only the distance per- 
pendicular to the nearest vortex filament is relevant, for 
a 3-dimensional quantum liquid we have n = 2, that is 
two degrees of freedom. The distance r of the test parti- 
cle from the vortex filament becomes a random variable 
given by 



(9) 



where again L is a suitable spatial scale introduced for 
dimensional reasons. 

We may estimate this length scale L as follows: For 
large distances r, of the order of average vortex filament 
distance d in the turbulent flow, the tracer particle fol- 
lows nearly normal type of Brownian motion, with Stokes 
law (|7|) valid in good approximation. Putting r = i into 
cq. (0), ® an d one arrives at the following estimate 
for the length scale L: 



L 
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(10) 



Clearly, our model requires small particles with a << d, 
if larger scales a > d are probed, one just get s ordinary 
Brownian motion with Gaussian behavior [28|, [29[ . 

The velocity distribution of the small tracer particle 
in the quantum turbulent flow described by eq.([2]) can 
now be calculated by using standard techniques of 
We first assume, for simplicity, a 
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which 



superstatistics 

constant 7 and define the parameter j3 
in equilibrium statistical mechanics corresponds to 
the inverse temperature, whereas here it is more a 
measure of distance from the nearest vortex filament. 
On time scales t satisfying 7 _1 <C t -C T a the stationary 
distribution of the tracer particle described by Eq.((2j 
for fixed (3 = =3r is given by the Gaussian distribution 



p{v\p) 



(11) 



assuming uniform distribution of the random vectors e. 
The situation becomes different for fluctuating /3, that is, 
if one allows the parameters 7 (or a) in Eq.@ to be vary- 
ing as well. Assuming that X\, . . . X n are independent 
Gaussian random variables, the resulting distribution of 
(3 = J27=i is a x 2 distribution of degree n, i.e. 
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The average of the fluctuating j3 is given by 

(/3) = n(X t ) = / 0f(fi) = p 
Jo 



(12) 



(13) 



and the variance by 



(14) 



The probability density to observe the velocity v of the 
test particle for any value of j3 is given by the marginal 
probability p(v) as follows 



P(v) = f(f3)p(v\(3)d(3 



(15) 



Substituting p(v\f3) and f(fi) from Eq.flTTJ and Eq.{l2 
into Eq. flS"]) . we obtain after a short calculation 



p(v) = 



r(f) 



(i + ^ 2 ) 



(16) 



These types of distributions play an important role in 
q-generalized versions of statistical mechanics [30], with 
the entropic index q related to the parameter n by q — 

n+l 

As we mentioned earlier the velocity of the tracer par- 
ticle depends on the perpendicular distance between the 
particle and the nearest evolving (and sometimes merg- 
ing) vortex filament. Therefore, the relevant degrees of 
freedom are n = 2 for 3-dimensional quantum turbulence. 
By substituting n = 2 in Eq. (fP6| one obtains 



p(v) = 



(2 + /?o« 2 )2 

Clearly, for large v this implies power-law tails 
p{v) cx v~ 3 . 



(17) 



(18) 



The remarkable result, however, is that we do not only 
get the power law tails but a concrete prediction for the 
entire shape of the probability distribution, including the 
region near the maximum. 

The probability distribution of kinetic energy E can be 
calculated from Eq. (116[) by using a simple transformation 
of random variables. For a particle of unit mass E = 
g(v) = i« 2 , hence v = g^ 1 {E) = \/2E and 



PE (E)dE = Pv (v) 2dv, 



(19) 



the factor 2 coming from the fact that there are two so- 
lutions ±v for the same energy E. This leads to the 
probability distribution of energy 



V2 



PE (E)=2 Pv (g-^E)) 

r(f) 



dg-\E) 



dE 



(20) 



E 



For n = 2 this predicts power law tails proportional to 
E~ 2 for large E. 
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FIG. 1: Experimental data of Paoletti et al [l[ and a fit using 
eq. (|2II) with variance parameter /9q = 4.5 and c = —0.12 for 
v x , respectively c = 0.54 for v z 



So far our model was based on a situation where the 
average velocity v of the particle is zero. Of course, in 
experiments there is often a drift velocity in the system 
that gives a non-zero mean velocity c to the test particle. 
In this case one has to replace v by v — c in the model 
equations we derived so far, and for n = 2 one ends up 
with 



p(v) 



{2 + fi {v-cfY 



(21) 



Let us now compare our model prediction with the ex- 
perimental data obtained by Paoletti et al. Fig. 1 
shows the experimental data for both velocity compo- 
nents v x and v z , and a fit by our analytic formula. An 
excellent fit is obtained. It is remarkable that the fit is 
not only correctly producing the power law tails but also 
the vicinity of the maximum. To illustrate this, Fig. 2 
shows the same data in a linear plot. 

Let us mention that our model directly predicts the 
power law exponent —3 in a universal way. The value —3 
is a consequence of the fact that vortex filaments are thin 
1-dimensional structures embedded into 3-dimensional 
space, thus leading to n = 2 in eq. ©. Our model also 
correctly reproduces the E~ 2 tails of the energy spectrum 
observed by Paoletti et al 0. 

Finally, we can also predicts the value (3q ~ 4.5 of 
the variance parameter to be used in eq. (|21[) . So see 
this, let us recall that Paoletti et al., in their experi- 
ment [l[, rescaled their measured velocity data to vari- 



ance 1. In these units their maximum velocity measured 
was v m ax — 10 (see Fig. 1). Strictly speaking, the vari- 
ance does not exist for any distribution that decays as 
v~ 3 for large \v\, but what exists is of course the vari- 
ance as calculated for a given experimental cutoff v max . 
From 




FIG. 2: Same as Fig. 1 but a linear scale is chosen, which 
emphasizes the vicinity of the maximum. The variance pa- 
rameter is slightly adjusted to obtain optimum coincidence in 
the central region. 



we obtain the predicted value /3q ~ 21ogu moa; s» 4.6, 
in agreement with what yields the optimum fit in Fig. 1. 
Thus, besides the (nonuniversal) systematic drift velocity 
c, all relevant parameters are predicted from first princi- 
ples. 

To conclude, in this paper we have developed a super- 
statistical dynamical model of Lagrangian quantum tur- 
bulence. This model predicts that the velocity statistics 
of small tracer particles in a quantum turbulent flow obey 
a power law distribution p(v) cx v ~ 3 and the distribution 
of energy follows a power law as well, i.e. p(E) oc E~ 2 . 
These results are in excellent agreement with Paolleti et 
al.'s measurements [l| as well as with the numerical re- 
sults obtained by other authors @, Q . Our theory pro- 
vides a universal prediction given by (f21~j) for both the 
center and the tail parts of the velocity distribution. The 
underlying stochastic model arises quite naturally out of 
the fact that small tracer particles see fluctuating effec- 
tive frictions, depending on the distance to the nearest 
vortex filament. 
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